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The apparent discrepancy between the Rosenbluth and the polarization transfer method for the 
ratio of the electric to magnetic proton form factors can be explained by a two-photon exchange 
correction which does not destroy the linearity of the Rosenbluth plot. Though intrinsically small, 
of the order of a few percent of the cross section, this correction is accidentally amplified in the case 
of the Rosenbluth method. 
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PACS numbers: 25.30.Bf, 13.40.Gp, 24.85. +p 

The electro-magnetic form factors are essential pieces 
of our knowlegde of the nucleon structure and this jus- 
tifies the efforts devoted to their experimental determi- 
nation. They are defined as the matrix elements of the 
electro-magnetic current J >1 {x) according to : 



< N{p')\J^{0)\N{p) > 
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where e ~ is the proton charge, M the nu- 

cleon mass, and Q 2 the squared momentum transfer. 
The magnetic form factor Gm is related to the Dirac 
{Fx) and Pauli (F 2 ) form factors by G M = Fi + F 2> 
and the electric form factor is given by Ge — F\ — ri*2, 
with r = Q 2 /4M 2 . For the proton, Fi(0) = 1, and 
i*2(0) = fi p — I = 1.79. In the one-photon exchange or 
Born approximation, elastic lepton-nucleon scattering : 



l{k)+N{p)^l{k')+N{p'), 



(2) 



gives direct access to the form factors in the spacelike 
region {Q 2 > 0), through its cross section : 



do B = C B {Q 2 ,£) 
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where e is the photon polarization parameter, and 
Cb{Q 2 ,s) is a phase space factor which is irrelevant in 
what follows. For a given value of Q 2 , Eq. © shows that 
it is sufficient to measure the cross section for two values 
of e to determine the form factors Gm and Ge- This is 
referred to as the Rosenbluth method 0]. The fact that 
da / C b{Q 2 , e) is a linear function of e (Rosenbluth plot 
criterion) is generally considered as a test of the validity 
of the Born approximation. 

Polarized lepton beams give another way to access the 
form factors |2(. In the Born approximation, the polar- 
ization of the recoiling proton along its motion {Pi) is 
proportional to G 2 M while the component perpendicular 
to the motion {P t ) is proportional to GeGm- We call 
this the polarization method for short. Because it is much 
easier to measure ratios of polarizations, it has been used 



mainly to determine the ratio Ge/Gm through a mea- 
surement of P t I Pi using [3j : 
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Thus, in the framework of the Born approximation, 
one has two independent measurements of R = Ge/Gm- 
In Fig. n we show the corresponding results, which we 
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which is common to both methods. The data are from 
Refs. 0,0,0- It i s seen that the deviation between the 
two methods starts around Q 2 = 2 GeV 2 and increases 
with Q 2 , reaching a factor 4 at Q 2 = 6 GeV 2 . A re- 
cent re-analysis of the SLAC cross sections [jj and new 
Rosenbluth measurements from JLab || confirm that 
the Rosenbluth and polarization extractions of the ratio 
Ge/Gm are incompatible at large Q 2 . This discrepancy 
is a serious problem as it generates confusion and doubt 
about the whole methodology of lepton scattering exper- 
iments. 

In this letter we take a first step to unravel this prob- 
lem by interpreting the discrepancy as a failure of the 
Born approximation which nevertheless does not destroy 
the linearity of the Rosenbluth plot. This means that we 
give up the beloved one-photon exchange concept and 
enter the not well paved path of multi-photon physics. 
By this we do not mean the effect of soft (real or virtual) 
photons, that is the radiative corrections. The effect of 
the latter is well under control because their dominant 
(infra-red) part can be factorized in the observables and 
therefore does not affect the ratio Ge/Gm- Here we con- 
sider genuine exchange of hard photons between the lep- 
ton and the hadron. Such higher-order corrections to the 
one-photon exchan ge a pproximation have been consid- 
ered in the past |9l llfj|. and their effects were found to 
be of order I - 2 % on the cross section. However, such 
estimates based on nucleon and resonance intermediate 
states can only be expected to give a realistic descrip- 
tion of the nucleon structure for momentum transfers up 
to Q 2 f GeV 2 , whereas they are largely unknown at 
higher values of Q 2 . 
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and their polynomial fits. 



Even if we restrict ourselves to the two-photon ex- 
change case, the evaluation of the box diagram (Fig. |2J) 
involves the full reponse of the nucleon to doubly virtual 
Compton scattering and we do not know how to perform 
this calculation in a model independent way. Therefore 
we adopt a modest strategy based on the phenomeno- 
logical consequences of using the full eN scattering am- 
plitude rather than its Born approximation. Though it 
cannot lead to a full answer it produces the following 
interesting results: 

• the two-photon exchange amplitude needed to ex- 
plain the discrepancy is actually of the expected 
order of magnitude, that is a few percent of the 
Born amplitude. 

• there may be a simple explanation of the fact that 
the Roscnbluth plot looks linear even though it is 
strongly affected by the two-photon exchange. 

• the polarization method result is little affected by 
the two-photon exchange, at least in the range of 
Q 2 which has been studied until now. 




FIG. 2: The box diagram. The filled blob represents the 
response of the nucleon to the scattering of the virtual photon. 

To proceed with the general analysis of elastic electron- 
nucleon scattering @, we adopt the usual definitions : 



k + k' 



and choose 

Q 2 = V, v = K.P, (6) 

as the independent invariants of the scattering. The po- 
larization parameter e of the virtual photon is related to 
the invariant v as (neglecting the electron mass m e ) : 



M 4 r(l 



M 4 r(l 
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For a theory which respects Lorentz, parity and 
charge conjugation invariance, the T-matrix for elas- 
tic scattering of two spin 1/2 particles can be ex- 
panded in terms of six independent Lorentz struc- 
tures which, following Ref. can be chosen 
as: u(k')u{k) u(p')u(p), u(k')u{k) u(p')j.Ku(p), 
u(k')^5u(k) u(p')75tt(p), u{k')^.Pu{k) u{p')^.Ku{p), 
u(k')j.Pu(k) u{p')u(jp), M(fc')757 .Pu{k) u{p')j^j.Ku(p). 
In the limit m e — > 0, the vector nature of the coupling 
in QED implies that any Feynman diagram is invariant 
under the chirality operation u(k) — > j^u(k),u(k') — > 
— u(fc')75. Therefore the Lorentz structures which change 
their sign under this operation must come with an ex- 
plicit factor m e . This allows us to neglect the structures 
which contain either u(k')u(k) or u(k')^u{k). Using the 
Dirac equation and elementary relations between Dirac 
matrices the linear combination of the remaining three 
amplitudes can be written in the form : 



x u{p') ( G M 7 
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where Gm, Fh, F 3 are complex functions of v and Q 2 , 
and where the factor e 2 /Q 2 has been introduced for con- 
venience. In the Born approximation, one obtains : 
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= G M (Q 2 ) 
= P 2 (Q 2 ), 
= 0. 



(9) 



Since P3 and the phases of Gm and F% vanish in the 
Born approximation, they must originate from processes 
involving at least the exchange of two photons. Relative 
to the factor e 2 introduced in Eq. (JHJ, we see that they 
are at least of order e 2 . This, of course, assumes that 
the phases of Gm and F 2 are defined, which amounts to 
supposing that, in the kinematical region of interest, the 
moduli of Gm and F2 do not vanish, which we take for 
granted in the following. Defining : 

G M = e^ M |G M |, F 2 = e^|F 2 |, F 3 = e*>\F 3 \, (10) 

and using standard techniques, we get the following ex- 
pressions for the observables of interest : 
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Pi 
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\G M \ ~ cos4>2m\F 2 \{1 



cos 4>?,m\Fz\p 



|Ga/|p + cos03m|^ 3 |(p 2 - r-r 2 ) 
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with 2 M = 4>2-4>M, 03Af = <t>3-4>M, 023 = 02 -03, and 

p = vjM 2 . If one substitutes the Born approximation 
values of the amplitudes I© then Eqs. (|11I12(I give back 
the familiar expressions of Eqs. 1)314(1 . 

To simplify the above general expressions, we make 
the very reasonable assumption that only the two-photon 
exchange needs to be considered. In practice we make an 
expansion in power of e 2 of Eqs. (|11I12|1 using the fact 
that 0M, 02 a- n d F 3 are at least of order e 2 but we do 
not expand \Gm\ and I-F2I, which is perfectly legitimate. 
This leads to the following approximate expressions : 



da 



C B (e,Q 2 ) 




(Ge/Gm) 2 with the coefficient of e, the two photon effect 
comes as a correction to a small number ~ (1/2. 79) 2 . So 
we expect that the correction will have a stronger effect 
in the Rosenbluth than in the polarization method. 

From Eqs. 1)131141) we see that the pair of observ- 
ables (da, Pt/Pi) depends on \Gm\, \Ge\, and TZ(F 3 ). 
In the first approximation, we know that \Gm(v,Q 2 )\ — 
G M (Q 2 ), \G e (v,Q 2 )\ ~ G E (Q 2 ), and only K(F 3 ) is re- 
ally a new unknown parameter. Thus allowing for two- 
photon exchange somewhat complicates the interpreta- 
tion of the lepton scattering experiments but not in a 
dramatic way. The main uncertainty is the dependence 
on v (or e) of F 3 and to further simplify the problem 
we make the following observations. First, if we look at 
the data of Ref. for da /Cb(£,Q 2 ) as a function of s 
we observe that for each value of Q 2 the set of points 
are pretty well aligned. We see in Eq. (|13f) that this 
can be understood if, at least in the first approximation, 
the product v F3 is independent of e. We do not have a 
first principle explanation for this but we feel allowed to 
take it as experimental evidence. To explain the linear- 
ity of the plot one must also suppose that \Gm\ and \Ge\ 
are independent of e (that is v ) but since the dominant 
term of these amplitudes depends only on Q 2 this is a 
very mild assumption. We then see from Eq. 1)13)1 that 
what is measured using the Rosenbluth method is : 
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where the neglected terms are of order e 4 w.r.t. the lead- 
ing one. By analogy, we have defined : 



G e = Gm-(1 + t)F 2 , 



(15) 



and 1Z denotes the real part. Note that Gf orn [y, Q 2 ) = 
Ge(Q 2 )- To set the scale for the size of the two-photon 
exchange term (F3) we introduce the dimensionless ratio : 



Y 2l {v,Q 2 ) =K 



vF, 



M 2 \G 



(16) 



M 



In the region of large Q 2 which is where the discrep- 
ancy really gets large, r is of order 1 or larger, while 
we can take as upper bound estimate |Ge|/|C?a/| ~ 
Ge(0)/Gm(0) = 1/2.79 . So, for a qualitative reason- 
ing, we can neglect \Ge\/\Gm\ with respect to r and, 
up to a term quadratic in Y27, the cross section has the 
form |G M | 2 (1 + e Y 2l ) 2 ■ So we expect Y 2l ~ a ~ 1/137. 
However in the Rosenbluth method where one identifies 



with \Ge\/\Gm\ and Y 2l essentially independent of e , 
rather than (R ,e ^ senbluth ) 2 = (Ge/Gm) 2 , as implied 
by one-photon exchange. Second, the experimental re- 
sults of the polarization method have been obtained for 
a rather narrow range of e, typically from e = 0.75 to 0.9 
for the points at large Q 2 . So, in practice, we can neglect 
the e dependence of Polarization and from Ec f C3} w e 
see that this experimental ratio must be interpreted as : 
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rather than Rp£ anzaUon - Ge/Gm- In order that 
Eq. 1)181) be consistent with our hypothesis we should find 
that Y 2l is small enough that the factor 2e/(l + e) intro- 
duces no noticeable £ dependence in R e poi ar i zat i n- 

We can now solve Eqs. 1)17118( 1 for |Ge|/|Gm| and Y 2l 
for each Q 2 . Since the system of equations is equiva- 
lent to a quadratic equation it is more efficient to solve it 
numerically. For this we have fitted the data by a polyno- 
mial in Q 2 as shown in Fig.^ and we shall consider this 
fit as the experimental values. In particular we do not 
attempt to represent the effect of the error bars which 
can be postponed to a more complete re-analysis of the 
data. The solution of Eqs. 1)17118)1 for the ratio Y^ p 
is shown in Fig. [21 where we can see that, as expected, 
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FIG. 3: The ratio Y£? p versus e for several values of Q 2 



it is essentially flat as a function of e and small, of the 
order of a few percent. Thus a tiny correction allows 
the Roscnbluth and the polarization method to give the 
same value for |Ge|/|GW |. It is reasonable to think that 
SGe = Ge—Ge and SGm — Gm~Gm are comparable to 
Y^ p and therefore |Ge|/|C?m| should not be very differ- 
ent from the actual value of Ge/Gm- So it makes sense 
to compare the value we get for -R l7 + 27 = |Gb|/|Gm| 

with the starting experimental ratios R^osenUuth an< ^ 
if' 1 1' 

Polarization ' 

that R?* +2l is close to R% p larization • 
tween the two curves can be attributed either to Y%* p or 



This is shown in Fig. 0] from which we see 
The difference be- 



to (SGm, SGe)- Insofar as (SGm, SGe) are of the same 
order of magnitude as Y^ p , which is small according to 
our analysis, our interpretation of this small difference 
is that the polarization method is little affected by the 
two-photon correction. 

In summary, the discrepancy between the Rosenbluth 
and the polarization method for Ge/Gm can be at- 
tributed to a failure of the one-photon approximation 
which is amplified at large Q 2 in the case of the Rosen- 
bluth method. The expression for the cross section also 
suggests that the two-photon effect does not destroy the 
linearity of the Rosenbluth plot provided the product 
TZ(vF 3 ) is independent of v. It remains to be investigated 
if there is a fundamental reason for this behavior or if it 
is fortuitous. Using the existing data we have extracted 
the essential piece of the puzzle, that is the ratio Y%* v 
which measures the relative size of the two-photon am- 
plitude F%. Within our approximation scheme, we find 
that Y^ p is of the order of a few percent. This is a very 
reassuring result since this is the order of magnitude ex- 
pected for two-photon corrections. What is needed next 
is a realistic evaluation of this particular amplitude. A 
first step in this direction was performed very recently in 
Ref. |l2j, where the contribution to the two-photon ex- 
change amplitude was calculated for a nucleon interme- 
diate state in Fig. [21 The calculation of Ref. 01 found 
that the two-photon exchange correction with intermedi- 
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with the value of fipR^Z^-y deduced from 



ate nucleon has the proper sign and magnitude to resolve 
a large part of the discrepancy between the two experi- 
mental techniques, confirming the finding of our general 
analysis. As a next step, an estimate of the inelastic part 
is needed to fully quantify the nucleon response in the 
two-photon exchange process. 

From our analysis we extract the ratio |Ge|/|Gm| 
which in the first approximation should not be very differ- 
ent from Ge/Gm- We find that it is close to the value ob- 
tained by the polarization method when one assumes the 
one-photon exchange approximation. This comparison 
is meaningful if, as suggested by the smallness of Y^ p , 
SGe and SGm are negligible. This could be checked by a 
realistic calculation of the two-photon corrections. How- 
ever we think that a definitive conclusion will wait for 
the determination of SGe and SGm as we did for Y^ p ■ 
The necessary experiments probably require the use of 
positrons as well as electron beams. 
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